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Abstract 

It  is  well  known  that  the  convex  hull  of  a  set  of  n  points  in  the  (Kudidcan)  plane  can  be 
found  by  an  algorithm  having  worst-case  complexity  ()(n  log  n).  In  this  note  we  give  a  short  linear 
algorithm  for  finding  the  convex  hull  in  the  case  that  the  (ordered)  set  of  points  from  the  vertices 
of  a  simple  (i.e.,  non-sclf-intcrsccting)  polygon. 
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1.  Introduction 

The  problem  of  finding  the  convex  hull  of  a  planar  set  of  points  P,  that  is,  finding  the 
smallest  convex  region  enclosing  P,  arises  frequently  in  computer  graphics.  For  example,  to  fit 
P  into  a  square  or  a  circle,  it  is  necessary  and  sufficient  that  H(P),  the  convex  hull  of  P,  fits;  and 
since  it  is  usually  the  case  that  H(P)  has  many  fewer  points  than  P  has,  it  is  a  simpler  object  to 
manipulate.  It  is  also  the  case  that  many  fast  graphics  algorithms  on  polygons  require  that  the 
input  polygon  be  convex,  thus  making  it  a  useful  preprocessing  step  sometimes  to  first 
transform  a  general  polygon  into  its  convex  hull.  A  number  of  algorithms  exist  for  finding  the 
convex  hull  of  a  set  of  points  (e.g.  [1  J[2][6]),  with  worst  case  complexity  0(n  log  n)  for  |P|  =  n.  It 
is  also  known  that  fl(n  log  n)  is  a  lower  bound  just  for  determining  H(P)  ••  that  is,  not  necessarily 
rendering  H(P)  in,  say,  clockwise  order  [7],  This  lower  bound  is  proved  for  a  decision  tree  model 
with  quadratic  tests,  which  accommodates  all  the  known  algorithms. 

An  interesting  case  of  the  convex  hull  problem  that  occurs  frequently  in  practice  is  when 
the  points  of  P  form  the  vertices  of  a  simple  polygon  (i.e.,  a  polygon  without  self-intersections). 
Several  authors  have  tried  to  find  a  fast  algorithm  for  this  problem. 'Sklansky  [5]  proposed  an 
O(n)  which  Bykat  [1]  later  showed  does  not  always  work.  It  has  also  been  noted  that 

the  algorithm  that  Shamos  [4]  suggested  can  sometimes  fail.  McCallum  and  Avis  [3]  recently 
published  an  O(n)  algorithm  which,  being  quite  complicated  and  utilizing  two  stacks,  entails 
rather  intricate  case  analysis  for  the  proof  of  its  validity.  In  this  note  we  provide  a  simple  linear 
time  algorithm  for  this  problem  and  a  short  proof  of  its  validity. 

2.  Preliminaries 

Let  P  be  a  simple  polygon  in  the  plane  with  n  vertices.  We  assume  that  each  vertex  vj  is 

represented  by  its  X  and  Y  coordinates.  Without  loss  of  generally  assume  that  P  =  <v^ ,  V£ . 

vn>  is  given  by  the  sequence  of  its  vertices  as  they  are  encou;  a  clockwise  traversal  of 

the  boundary.  (The  orientation  of  P  can  be  easily  tested  and,  if  n<*c.  ry,  reversed.)  The 
convex  polygon,  also  clockwise  oriented,  whose  vertices  are  the  set  of  extreme  points  of  P  will 
be  denoted  by  H(P).  For  a  polygon  P,  we  will  use  P[vj,  Vj]  to  denote  the  path  in  P  from  Vj  to  vj 
(following  the  orientation  of  P).  The  concatenation  of  two  paths  p  and  q  will  be  written  as  p«q. 
Given  two  points  x  and  y,  the  notation  L[x,  y]  will  refer  to  the  (directed)  line  segment  from  x  to  y. 
The  primitive  test  used  in  our  convex  hull  algorithm  is  a  right  turn  test.  A  turn  is  an  ordered  triple 
of  three  points  (x,y,z).  We  say  that  (x,y,z)  is  a  right,  straight  or  left  turn,  written  (x.y.z)  =  1,0  or 
- 1,  if  the  Z  component  of  the  cross  product  (z-  y)X(y-  x)  is  positive,  zero  or  negative, 
respectively.  (See  Figure  1 ).  Alternatively,  we  say  that  z  lies  to  the  right  of,  on  or  to  the  left  of 
L[x,y]  if  (x,y,z)  is  a  right,  straight  or  left  turn,  respectively.  (One  can  test  whether  P  is  clockwise 
oriented  by  finding  a  vertex  Vj  with  minimum  X  coordinate  and  checking  that  (vj.-j ,  vj,  vj  +  -j)  is 
not  a  left  turn.) 


Let  the  vertices  of  the  convex  hull  H(P)  be  <q^ , ....  q^>.  It  follows  from  the  Jordan  Curve 
Theorem  that  <q-| , qh>  corresponds  to  a  subsequence  of  the  vertices  in  P,  with  the  following 
characterization:  (Figure  2) 

(A1)  <q-| . qii>  forms  a  convex  polygon  (clockwise  oriented). 

(A2)  No  point  of  the  path  P[qj.i,qj]  lies  to  the  left  of  l[qj.-|,qj]  for  2<i<h+1.  (Hereqh  +  iis 

-identified  with  q-|). 

We  will  call  a  polygon  L[qj.i,  qj]oP[qj,  qj.-j]  which  satisfies  (A2)  a  pocket,  denoted  by  Kj,  and 
edge  Lfaj.-j ,  q,]  the  top  of  the  pocket.  Notice  that  points  of  Pfoj,  qj.i],  other  than  its  two  end 
points  qj  and  qj.-j ,  may  lie  on  L[qj.-j ,  qj].  Thus,  in  general,  each  pocket  Kj  has  an  interior  which  is 
a  disjoint  union  of  a  finite  number  of  regions  lying  to  the  right  of  Lfaj.-j ,  qj];  the  interior  may  also 
be  empty  in  the  degenerate  case  when  P[qj.-j ,  qj]  *  L(q,.i ,  qj].  The  task  of  finding  H(P)  can  then 

be  specified  as  finding  a  sequence  of  pockets  R^ ,  R2 . Rh  so  that  their  tops  L[q-j ,  q2],  L[q2. 

q3], ....  Lfah,  q-j]  form  a  convex  polygon.  Now,  suppose  we  specify  two  vertices  of  H(P)  --  say, 
q-l  and  qs.  Then  the  chord  L(q1 ,  qs]  divides  H(P)  into  two  convex  polygons  with  L[q-| ,  qs]  as  a 
common  edge,  while  their  other  edges  satisfy  (A2).  Indeed,  (At)  and  (A2)  above  imply  that  a 
characterization  H(P)[q-|,  qs]  =  <q-|, ....  qs>  is  as  follows:  (Figure  2) 

(B1)  <q-| , qs>  forms  a  convex  polygon  (clockwise  oriented). 

(B2)  L[qj_i,  qj]<>P[qj,  q,.-|]  forms  a  pocket  for  2<i<s. 

Note  that  H(P)[q-|,  qs]  depends  only  on  P[q-|,  qs]-  We  will  call  HfP)^,  qs]  the  left  hull  of 
P[q1 ,  qg].  In  practice,  a  convenient  choice  for  q-|  and  qs  could  be,  say,  two  points  of  P  with  the 
minimum  and  the  maximum  Xcoordinate,  respectively.  Without  loss  of  generality,  we  now 

assume  that  in  P  =  <v1 ,  v2 . vm,  vm  +  . vn>,  vertices  v-j  and  vm  have  the  minimum  and  the 

maximum  X- coordinate,  respectively. 

The  problem  of  finding  H(P)  thus  reduces  to  that  of  finding  the  left  hulls  of  P]v^ ,  vm]  and 
p(vm,  v-)].  We  will  describe  an  algorithm  for  the  latter  problem  in  the  next  section. 

3.  The  Algorithm 

Algorithm  LeftHull  is  given  by  the  diagram  in  Figure  3.  The  algorithm  considers  the  points 
<v-j ,  v2, ....  vm>  of  P[vi,vm]  in  order.  The  main  data  structure  used  is  a  stack  Q  =  <q$,  qi, .... 
qt>,  with  q$  being  the  element  at  the  bottom,  and  the  variable  t  pointing  to  the  top  of  the  stack. 
(The  element  q$  corresponds, to  qs  of  the  preceding  discussion.)  The  main  body  of  the 
algorithm  is  a  loop,  consisting  of  two  boxes  N  and  C. 
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Algorithm  LcftHull 
Input: 

P(v1  •  vml  =  <V1 vm> 

Output: 

HCPJIvt  ,  vm] »  Q  =  <q$,  qi . qt> 


Figure  3.  Algorithm  LeltHull. 


Each  time  we  enter  box  N,  the  following  induction  hypotheses  are  true  (see  Figure  4). 

(i)  Q  =  <q$,qi . qt>  forms  a  convex  polygon  (clockwise  oriented),  with  q$  =  vm,  and 

qi=vi. 

(ii)  L[qj.i ,  qj]oP[qj,  qj.-|]  forms  a  pocket-for  2<i<t. 

The  only  difference  between  these  induction  hypotheses  and  the  characterization  of 
H(P)[qi,  q$l  given  by  (B1)  and  (B2)  is  that,  here,  no  stipulation  is  made  about  L[qt,  q$]  and  Pfat. 
q$]  in  statement  (ii).  Of  course,  we  must  show  that  L[qt,  q$]oP[q$,  qt]  indeed  forms  a  pocket 
when  Algorithm  LeftHull  terminates.  To  explain  the  algorithm,  we  first  define  the  notation  used. 

1 .  x  refers  to  the  input  point  Vj  under  consideration,  and  y  refers  to  the  input  point  vj 
preceding  the  point  q^ 

2.  We  use  b  to  denote  the  turn  (y,qt  x);  and  c  to  denote  the  turn  (qt.-|  qt,x).  For  the  current 
input  x  to  be  added  to  the  stack,  condition  c  =  1  is  required  by  (i).  The  function  of  b  will 
be  explained  later. 

3.  "Pushing  x"  means  executing  [t«-t+  1;  qt*-x;  update  y];  "popping  x"  means  [t«-t-1].  No 
code  is  executed  for  "rejecting  x". 

4.  By  "x«-  input"  we  mean  [IF  input  is  not  exhausted  THEN  x  *■  nextPoint  ELSE  EXIT]. 

We  will  go  through  Algorithm  LeftHull  once,  explaining  the  major  steps  and  proving 
induction  hypotheses  (i)  and  (ii).  In  box  I,  we  do  the  initialization  and  form  K-j  by  picking  the  first 
point  after  q^  with  c  =  1  to  be  Q2-  Thereafter,  each  iteration  of  boxes  N  and  C  creates  a  new 
pocket  while  maintaining  overall  convexity.  Inductively,  when  we  enter  box  N,  pockets  K-j, K( 
satisfying  (i)  and  (ii)  have  been  found,  with  x  being  the  first  input  point  seen  after  qt.  We  will 
consider  two  cases,  corresponding  to  whether  b  =  (y,qtix)>0  or  not. 

Case  (a).  Suppose  b>0  (Figure  5a).  The  DO  UNTIL  loop  in  box  N  finds  the  first  point  x  to  the 
left  of  L[q{,  q$].  This  defines  a  new  pocket  bounded  by  L[qt,  x)oP[x,  qj. 

Case  (b).  Suppose  b  =  -1  (Figure  5b).  Then  x  could  be  either  in  the  pocket  Kj,  or  outside  of  it. 

Note  that  in  the  latter  case  P  must  be  left  of  L[qt.i ,  q^.  In  the  former  case,  P  must  also 
eventually  leave  Kt  by  moving  left  of  L[qt_i ,  qt]  at  some  point  x  ••  since  this  is  the  only 
way  the  simple  polygon  P  could  reach  q$  without  crossing  P[qt.i ,  qj.  When  such  a 
point  x  is  found,  we  will  make  a  new  pocket  out  of  by  making  L[q^i ,  x]  its  top.  This 
is  accomplished  by  the  last  statement  in  box  N  which  pops  qt  from  the  stack.  (In  fact, 
this  statement  is  included  in  box  N  purely  for  ease  of  discussion;  it  is  clearly 
superfluous  in  view  of  the  first  statement  of  box  C.) 


Thus,  in  either  case  (a)  or  (b),  after  executing  box  N  we  will  have  t  +  1  pockets  lying 
between  q-j  and  x.  They  would  satisfy  the  induction  hypotheses  if  c  =  1  were  true.  Box  C  is  then 
executed  to  fulfill  this  convexity  requirement,  by  merging  the  last  pocket  with  as  many  pockets 
as  necessary  from  the  stack  until  the  condition  c  =  1  is  satisfied.  (This  will  happen  before  q-j  is 
popped,  since  (q$,  q-| ,  x)  is  a  right  turn  as  one  can  easily  show.)  At  this  time,  x  becomes  the  top 
element  of  Q,  and  the  inductive  step  is  complete. 

The  preceding  discussion  can  easily  be  made  rigorous.  The  correctness  of  the  algorithm 
thus  only  depends  on  extending  assertion  (ii)  to  include  the  last  segment  P[qt,  q$]  when  the 
algorithm  terminates.  If  qt  itself  is  the  last  input,  then  (ii)  is  trivially  true.  Otherwise,  the 
algorithm  must  terminate  while  executing  the  DO  statement  in  box  I  or  R.  As  we  argued  in  the 
preceding  paragraph,  it  is  not  possible  that  the  algorithm  terminates  in  box  R  after  having 
entered  it  with  b  =  0.  Thus  it  must  terminate  either  in  box  I,  or  in  box  R  after  entering  it  with  b  =  1. 
In  either  case,  the  property  "x  lies  to  the  right  of  L[q{,  q$]"  will  be  true  for  all  input  points  x  seen 
after  q^.  This  completes  the  argument.  The  algorithm  is  obviously  linear,  since  each  input  point 
can  be  pushed  onto  or  popped  from  the  stack  at  most  once,  if  it  is  not  rejected  outright. 

4.  Conclusion 

We  present  a  simple  linear  time  algorithm  for  finding  the  convex  hull  of  a  simple  polygon. 
Note  that  our  algorithm  can  actually  be  applied  to  non-simple  polygons  P  as  follows:  First,  at 
each  point  p  of  intersection  of  two  (or  more)  edges  of  P,  place  a  new  vertex  v(p).  Viewed  as  a 
graph  G(P),  all  vertices  (old  and  new)  have  even  degree,  so  that  the  graph  G(P)  is  Eulerian. 
Choose  a  fixed  Eulerian  circuit  in  G(P).  It  is  not  difficult  to  see  that  each  new  vertex  v(p)  of 
degree  2d(p)  may  now  be  split  into  d(p)  slightly  perturbed  vertices  Vj(P),  1  <i<d(P),  of  degree  2 
so  that  each  vj(P)  is  in  the  interior  of  H(P)  and  the  resulting  polygon  P*  ( =  Eulerian  circuit)  is 
simple.  Therefore  H(P*)  =  H(P)  and  our  algorithm  can  be  used  to  compute  H(P*).  Observe, 
however,  that  if  P  has  n  vertices,  P*  can  have  cn^  vertices. 

In  the  case  that  P  is  a  simple  path  with  n  vertices,  we  can  join  the  first  and  last  points  of 
P  creating  a  simple  new  edge  to  form  a  (possibly  non-simple)  polygon  Pi  By  applying  the 
preceding  transformation  we  can  then  find  a  simple  polygon  P*  on  at  most  3n  vertices  with 
H(P*)  =  H(P)  =  H(P)  in  time  0(n). 
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